A Spallation Model for the Titanium-rich Supernova Remnant Cassiopeia A 
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Titanium-rich subluminous supernovae are rare and challenge current SN nucleosynthesis models. 
We present a model in which ejecta from a standard Supernova is impacted by a second explosion 
of the neutron star (a Quark-nova), resulting in spallation reactions that lead to 56 Ni destruction 
and 44 Ti creation under the right conditions. Basic calculations of the spallation products shows 
that a delay between the two explosions of ~ 5 days reproduces the observed abundance of 44 Ti in 
Cas A and explains its low luminosity as a result of the destruction of 56 Ni. Our results could have 
important implications for lightcurves of subluminous as well as superluminous supernovae. 

PACS numbers: 23.23.+X, 56.65.Dy 



I. INTRODUCTION 

Cas A is a young galactic supernova remnant formed 
in the aftermath of a core-collapse explosion of a mas- 
sive star. Tentative classification as Typellb (based on 
the detection of weak Helium lines in addition to the 
Hydrogen line) has been strengthened by infra-red stud- 
ies of the scattered light echo [l|. Cas A is intensively 
studied as a prototype that can unveil the pristine com- 
position of nucleosynthetic yields, thereby constraining 
aspects of the supernova mechanism and stellar evolu- 
tion models. Multi-wavelength studies show that emis- 
sion from the remnant is dominated by a bright ring, but 
also find jets and knots, as well as an X-ray bright CCO 
(central compact object) believed to be a rapidly cool- 
ing neutron star. Cas A is an unusual supernova (SN) 
in some respects: nuclear decay lines of 44 Ca (1157 keV) 
and 44 Sc (67.9 and 78.4 keV) detected by COMPTEL 
and BEPPO SAX indicate a very large synthesized 44 Ti 
mass of (0.8-2.5) xl(T 4 M (M Q is the solar mass), 
which would imply an ejected 56 Ni mass of at least 0.05 
Mq 0, 0| , making Cas A extremely bright given its prox- 
imity; but no definitive historical record of such a bright 
SN at that time (~ 300 yrs ago) exists. The possible 
detection by Flamsteed in 1680 suggests 6th magnitude. 
Possible solutions center on extinction due to a surround- 
ing dust cloud, possibly generated by a pre-SN wind from 
the massive star and an asymmetric SN explosion Q . Al- 
though observations of Cas A do give indications of these 
features, they do not explain why no other Ti-rich SNc 
have been found despite searches in massive star regions 
in the inner galaxy @, 0| • It appears that Ti-producing 
SNe are quite rare. NuSTAR (nuclear spectroscopic tele- 
scope array) aims to map the 44 Ti in more Cas A-likc 
remnants to solve this puzzle. 

In this letter, we present an alternative that can rec- 
oncile the sub-luminous nature of Cas A with its excess 
44 Ti production. Assuming that the progenitor of the 
Cas A remnant was not atypical (i.e., did not produce 
large amount of Ti in situ) , we suggest that Ti is formed 



instead as a spallation product when neutron-rich mate- 
rial ejected from a "second" explosion , viz., that of the 
neutron star (a quark- nova), impacts and destroys the 
56 Ni layer ejected by the preceding explosion (the SN). 
The basic picture is that a SN can produce a massive neu- 
tron star, which then converts explosively to a quark star 
(QS hereafter) (in an event called a Quark Nova or QN; 
[SHI)- Such an explosion can happen if the Neutron Star 
(NS), in its spin-down evolution, reaches the quark de- 
confinement density and subsequently undergoes a phase 
transition to the more stable strange quark matter phase 
[T0I - IT2I . resulting in a conversion front that propagates 
toward the surface in the detonative regime. The outer 
layers of the parent NS are ejected from an expanding 
thermal fireball [IH, [l4| which allows for ejecta with ki- 
netic energy easily exceeding 10 52 erg. If the QN occurs 
less than a few weeks after the SN, the QN energy re- 
lease reheats the preceding SN ejecta, and results in an 
event we call a dual-shock Quark Nova (dsQN hereafter). 
In previous papers, we introduced the dsQN as a model 
for superluminous SNe [15| and discussed their photo- 
metric and spectroscopic |16| signatures. Here, for the 
first time, we explore the nuclear processing of the inner- 
most SN ejecta by the QN relativistic ejecta (neutrons 
and heavy nuclei). We show that such a model for Cas 
A can explain why the SN is sub-luminous yet Ti-rich; 
in addition, the rarity of such events follows naturally 
as a constraint from the time delay between the two ex- 
plosions. Testable predictions based on our model are, 
delayed Hydrogen signatures (weeks after the second ex- 
plosion) and a modified SN light curve. 



II. THE SPALLATION MODEL 

Beam and Target: In analogy with spallation reactions 
in the laboratory, we frame our model in the context of 
a "beam" and a "target". The QN provides the "beam": 
a relativistic outflow of neutron-rich material from the 
NS surface, caused by an explosive phase transition in 
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its core to a more compact quark phase. Recent nu- 
merical studies of the phase conversion front [TtJ sug- 
gest supersonic laminar motion of the conversion front, 
which can become unstable [l8| , wrinkling the conversion 
front to serve as a platform for a DDT (deflagaration-to- 
detonation). The outcome, depending on the conversion 
efficiency of the shock to kinetic energy, is ejection of 
about M QN ~ 10 _3 M Q HQ of the NS's outermost layers 
at nearly relativistic speeds with average Lorentz factors 
ofTQN ~ 10. The total number of ejected nucleons in this 
beam (mostly neutrons) is then TV ~ 1.2 x 10 54 Mqn _ 3 
where Mqn.-3 is the QN ejecta mass in units of 1O -3 M0. 
Adopting these fiducial values, the neutron energy is 
Eq ~ 10 GcV. This beam of relativistic neutrons (speeds 
close to c), will overtake and strike the innermost layers 
of expanding SN ejecta (the "target") which are mov- 
ing much slower at a speed " v <C c" . Setting our clock 
by the SN explosion at t=0, this collision will happen a 
time ideiayj the delay between the SN and the QN explo- 
sions. The collision between the QN and SN ejecta causes 
spallation and subsequently other nuclear reactions. The 
crux of our argument is that these spallation reactions 
can be a mechanism to explain some of the unique fea- 
tures of Cas A discussed previously, if idciay is chosen 
appropriately. 

We assume an onion-like profile of the expanding 
shocked SN ejecta (i.e., no mixing) with the innermost 
ejecta, viz., Ni nuclei (mass number A=56) constituting 
the target at a distance from the CCO of Rj n (t)—v ideiay 
The target number density in the Ni layer is approxi- 
mately constant at ua=Ma/ (4:irRf n AR), where AR is 
the thickness of the Ni layer. The neutron mean free path 
for spallation in the Ni layer is A = 1/ (riACsp), where the 
spallation cross-section for neutrons on a target nucleus 
is empirically described with a sp ~ 45A - 7 f(A) mb with 
f(A) ~ 1 a factor of order unity (l9^ . We have ignored a 
weak dependence on energy (< 10% for energies E > 100 
MeV; 19]), so that spallation mean free paths for the 
few subsequent generations in the spallation cascade are 
roughly constant. The average number of collisions an 
incoming neutron makes in the Ni layer is 
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where Ma is in units of 0.1M Q . The SN ejecta must be 
sufficiently dense that N co \\, > 1, which for fiducial values 
of Ma and v in Eq.([T]), limits idciay < 8.3 days. 

Spallation Reactions: We include only the two most 
relevant reactions: spallation by neutrons and by protons 
(e.g. n + A — > products + n x (£ nn + 1) + p x £ np ) where 
£ is the multiplicity. The n + A neutron multiplicity as 
a function of the beam energy and target material shows 
roughly linear dependence on the target mass number (in 
the range 12 < A < 238). We use the semi-empirical for- 
mula ( nn {E) ~ A (0.0833 + 0.03171ni?), where the neu- 
tron energy E is in GeV[20]. This formula gives better 
than 10% accuracy for A > 40. The average total multi- 



plicity is 

C(E,A) = Cn„ + Cn P - 4.67A 56 (1 + 0.381n£)r np , (2) 

where Y np = (1 + Cn P /Cnn) is in the range 1.25 < Y np < 
1.67 (e.g. [20]). For proton induced reactions (p + A), 
most of the results are similar to those of neutron induced 
reactions [20[ • We treat neutrons and protons identically. 
Spallation effectively ceases when ( drops below 1, cor- 
responding to projectile energy of E ~ 73 MeV. 

Spallation Statistics: Let us divide the Ni layer into 
~ iV co ii imaginary sub-layers of radial thickness ~ A, with 
k=0 denoting the innermost one. A given Ni nucleus in 
this layer will be hit multiple times (-/Vhits) by neutrons, 
resulting in a product nucleus 
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E° is the typical energy of neutrons impinging on the 
k=0 layer. To produce a realistic distribution of product 
nuclei from this sub-layer, we draw N® its and £° from 
Poisson distributions peaking at A^°hits and £°, where 
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and (°=((E° , Ai). For subsequent sub-layers, it follows 
that 
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where r\ accounts roughly for the incident energy removed 
by radiation, nuclear excitation and other sub-products 
(pions, He, Deuterium etc ...). Since t he p ions and other 
multiplicities are small, we set r\ = [20t j . A Poissonian 
description is appropriate since we have a small spalla- 
tion cross-section but a rapidly increasing -/Vhits due to 
the cascading effect. The total number of layers that 
experience spallation is given by min (-/V co n, fc max ) where 
fi 1 """ 1 ~ 73 MeV, since for a specified thickness (mass) 
of target material, one can run out of material before 
spallation becomes insignificant. 



III. RESULTS 

Spallation products in the inner SN ejecta: Fig. 1 
shows the probability distribution of the spallation prod- 
ucts from 56 Ni layer for fdeiay of 4, 5 and 6 days. 56 Ni is 
depleted, while 44 Ti and light elements (H through Ne) 
are produced. The ideiay = 5 day Ni target case is partic- 
ularly interesting for Cas A: note the Ti and C production 
peaks. Figure 2 shows normalized mass yields of spal- 
lated fragments r/A = Ma/Ma , where Ma is the initial 
amount of Ni. We observe that for 3 days < idoiay < 7 
days, the result is a Ni-poor, Ti-rich, C-rich debris. 



3 



(A = 56, < delay = 4 days) (A = 56, t ielay = 5 days) (A = 56, < delay = 6 days) 
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FIG. 1: Spallation products in successive layers (back to front) from 56 Ni for idciay =4,5,6 days. The overall distribution is 
shown in the front layer labelled "combined". 



Neutron- capture versus neutron- decay: What happens 
to the neutrons themselves? For E k " m < 73 MeV fol- 
lowing the spallation regime, the neutrons mainly de- 
posit heat through inelastic collisions. Neutron-capture 
happens once the neutron energy is further reduced to 
< 30 MeV. Compared to the free neutron lifetime of 
~ 720 seconds, the rather long n-capture timescale of 
T cap. = l/( ri A <T cap. w n) ~ 10-2 hours (for fiducial values of 
v u ,Ma, AR where v n is the neutron thermal speed and 
Ocap. the capture cross-section [21| ) implies that ideiay 
must be shorter than about 16 hours for significant neu- 
tron capture. In such cases, nucleosynthetic yields can 
be slightly altered from those expected in a SN. 

Hydrogen Formation: There are three sources of H 
in our model: (i) a direct byproduct of spallation (i.e. 
when A is reduced to 1) for short time delays, which 
contributes at most 1% of the total mass in the target 
layer; (ii) from spallation protons forming hydrogen via 
recombination once their energy E ~ 73 MeV (protons 
are about 1/3 of the total nucleons formed); (hi) from 
proton recombinations following /3-decay of the neutrons 
that evade capture: for idciay > 0.7 days, most neutrons 
decay to protons which form H by recombination. The 
total amount of hydrogen that we estimate could form 
from these three factors is Mh ~ 0.1M©Mqn,-3- 



IV. CONCLUSIONS AND PREDICTIONS 
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FIG. 2: Mass Yields (upper panel) and total spallation neu- 
trons&protons (lower panel) versus time delay for 56 Ni target. 
In the lower panel, spallation layers are numbered 0-9, with 
spallation effectively ceasing after 10 layers for tdoiay = 2 days. 



Ni-poor, Ti-rich SNe: We suggest that the Ti-rich Ni- 
poor yield of Cas A can be explained if the SN was fol- 
lowed by a QN with ideiay ~ 5 days. From Fig.s 1&2, 
this reproduces the observed abundance trends of key ele- 
ments. In our model, 44 Ti is formed from the destruction 
of 56 Ni and thus Ti-rich dsQNe will be necessarily sub- 
luminous. Also, the rarity of Ti-rich SNe in the massive 
star populations could be because the mass-cut is above 
the Ti zone in these SNe, suppressing Ti ejection. Only 
those experiencing a QN explosion following the SN, with 



the appropriate tdciay, would show Ti produced by spal- 
lation (as we suggest for Cas A). The higher mass-cut 
would also favor more massive NSs, ideal candidates for 
the QN transition to occur 22] . 

Delayed Hydrogen signatures: Decay of spallation neu- 
trons to protons will not immediately form hydrogen 
since the recombination and thermal continuum radia- 
tion remains trapped and ionizes the H until it escapes 
when Ai?A = l/( n c,Ao°Th.) where <7xh. is the Thompson 
cross-section and n o .A the electron density in the target 
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layer. For Aq = 56, the recombination radiation escapes 
when the innermost layers of the SN ejecta reaches a ra- 
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hydrogen should be observed ~ 41 daysX M A ' j/vin.sooo 
days following the explosion. 

Other signatures: The late time SN lightcurve should 
carry signatures of 56 Ni and 56 Co destruction by the QN 
ejecta. In the case of 56 Ni an amount of it would have 
decayed in idciay and converted to 56 Co, so that, for eg., 
only 57% of the original 56 Ni will experience spallation 
from the QN for idciay = 5 days. The remaining 43% of 
the original 56 Ni would have decayed contributing mostly 
to PdV work in the SN ejecta. This means that the effi- 
ciency of the QN spallation determines how much of the 
remaining 56 Ni is destroyed and directly influences the 
SN luminosity. In general, we find that for idciay < 8 
days, dsQNe appear as Ni-poor (sub-luminous) and Ti- 
rich/C-rich SNe. If idciay ~ 8 days, the Ni decay lumi- 
nosity should be only slightly reduced but the Co decay 
luminosity will be strongly suppressed (i.e. no 77 day 
tail). Even smaller idciay can result in earlier destruc- 
tion of 56 Ni and thus lead to a subluminous SN with 
no 77 day Cobalt tail. It is interes ting to consider the 
case of a Si layer (with ~ 10% Ca; |23|) over a Ni layer 
where the neutrons exit the Ni layer with E n > > 73 MeV 
(i.e. fc max >> iVcoii). For fiducial values, this occurs for 
^dciay ~ 6-8 days (iV co u ~ 2; see Fig. 2, lower panel). 
This implies that the Si/Ca is destroyed by spallation, 
which could explain the lack of lines such as Si, Ca and 
Iron in some superluminous SNe j24l] . This delay is also 
ideal for explaining reenergization of superluminous SNe 
15]. A <dciay ~ 5 day dsQN leads to C formation (see 
Fig. 2); this may explain the C-rich atmosphere of the 



compact object in Cas A [25J, which, in our model, would 
be a radio-quiet QS (an aligned rotator; (2(|) surrounded 
by a layer of fallback material rich in C. Finally, the QN- 
ejecta also contains heavy elements from neutron-capture 
(see [27]]) and collisions between these and the inner lay- 
ers of the SN ejecta will lead to other unique signatures 
which we have yet to investigate. 

Model assumptions: Our model has some fine tuning 
that is unavoidable due to uncertainties regarding the 
nature of the hadron-quark phase transition. 25[ give 
mass-radius constraints for the compact object in Cas A, 
which effectively rule out low-mass QSs based on non- 
interacting quark equations of state. However, large and 
heavy QSs may exist, so long as the quark superconduct- 
ing gap and strong coupling corrections are taken into 
account (e.g. |28l - [30l ] ) . The issue of the mass-radius re- 
lation for quark stars is still a matter of debate. Spectral 
fitting of Cas A agrees very well with theoretical cooling 
models for NSs, when superfluidity and pair-breaking ef- 
fects are taken into account [3l|. It is unlikely that a 
QS would exhibit exactly the same cooling behaviour as 
a NS, which is a problem for our model, but there is 
no comprehensive cooling simulation studies of QS and 
it might be purely coincidental that at this particular 
young age, a NS and a QS have the same surface tem- 
perature. Studies of cooling of QSs which include similar 
attention to physics details (e.g. color superconductiv- 
ity) are needed to determine whether they could be at 
all consistent with Cas A. 
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